i 


iKSH 


'■■  111  1 


■■■■■.•■'■■ 


KfflhvD 


.;,.  ; 


-H^M 

,;; f'7;,:\;,L ■ .  ■ 

.Maw 


Dudltv      .ibrary,  NPS 
Monterey,  CA  93943 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


A  LINEAR  APPROXIMATION  OF  THE  SOURCE 
POSITION  USING  MULTIPLE  MAD 

by 

Wolf -Hub ertus  Bock 


September  1983 


Thesis  Advisor:   Andrew  R.  Ochadlick,  jr 


Approved  for  public  release;  distribution  unlimited 


T210144 


UNCLASSIFIED 


Dudley  Knox  Library,  I 

Wnn+  orav 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Bntarad) 


REPORT  DOCUMENTATION  PAGE 


r     REPORT   NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO 


3.     RECIPIENT'S  CATALOG  NUMBER 


4.    TITLE  fand  Subtltla) 

A  Linear  Approximation  of  the  Source 
Position  Using  Multiple  MAD 


5.     TYPE  OF   REPORT  a   PERIOD  COVEREO 


Master's    Thesis 
Sppt-Pmhpr    1983 


£ 


S.     PERFORMING  ORG.    REPORT  NUMBER 


7.    AUTHOR<a> 

Wolf-Hubertus   Bock 


8.     CONTRACT  OR  GRANT  NUMBERC*.) 


•      PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


tt.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93943 

14.    MONITORING  AGENCY  NAME  ft   AOORESSfU  dlllarant  from  Controlling  Otflca) 


12.     REPORT   DATE 


September    1983 


13.     NUMBER  OF   PAGES 


61 


15.     SECURITY  CLASS,   (ol  thta  raport) 


15a.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


IS.     DISTRIBUTION  STATEMENT  (ol  thia  Raport) 

Approved  for  public  release;  distribution  unlimited 


17.    DISTRIBUTION  STATEMENT  (ol  tha  abstract  antatad  In  Block  30,  It  dlllarant  from  Raport) 


IS.    SUPPLEMENTARY  NOTES 


IS.     KEY  WORDS  (Cotulnua  on  ravaraa  aid*  II  nacaaaary  and  Idantlty  by  block  numbar) 

MAD  localization,  Multiple  MAD,  Non-acoustic  localization 


20.     ABSTRACT  (Coitlirtua  an  ravr—  aid*  II  nacaaaary  and  Idantlty  by  block  numbar) 

For  certain  assumptions,  an  analysis  of  multiple  MAD  signals 
results  in  a  reasonable  estimate  for  the  localization  of  a 
target  relative  to  the  MAD  platform.  This  is  achieved  by  using 
selective  approximations  to  linearize  an  initially  nonlinear 
problem.  The  simulation  ignores  noise  and  requires  an  estimate 
of  the  magnitude  of  the  target  magnetic  moment  components . 
Results  indicate  that  the  best  localization  estimates  are 


0D 


I  JAN  7J     W J 


EDITION  OF   1  NOV  «S  IS  OBSOLETE 
S/N  0102-  LF-0U-  6601 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF   THIS  PAGE  (Whan  Data  Bntarar 


u  Li  \j  1-tJn.Lj  i~j  j-j.    j.i^±s 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dmtm  Bnfnd) 


achieved  when  the  platform  is  on  cardinal  headings ,  and  when 
the  target  moment  has  a  strong  vertical  component. 


S-N  0102-  LF-  014-  6601 


JNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEr*h»n  Dttm  Bnfrmd) 


Approved  for  public  raiease;  distribution  unlimited 


&  Linear  Approximation  of  the  Source 

Position  Using  Multiple  HAD 


by 


Wolf-Huber-us  ^ock 
Lieutenant,    United   Statas   Navy 
3.  A.,   Rice    University,    1975 


Submitted    in  partial   fulfillment    of   the 
requirements    for  the    dagree   of 


MASTER    3?    SCIENCE    IN    SYSTEMS    TECHNOLOGY 
(ANTISUBMARINE    WARFARE) 


from   the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1933 


Mo;  $5%8fS 


: 


ABSTRACT 


Fcr  certain  assumptions,  an  analysis  of  multiple  MAD 
signals  results  in  a  reasonable  estimate  for  the  localiza- 
tion of  a  target  relative  tc  the  MAD  platform.  This  is 
achieved  by  using  selective  approximations  to  linearize  an 
initially  nonlinear  problem.  The  simulation  ignores  noise 
and  requires  an  estimate  of  the  nagnitude  of  the  target 
magnetic  moment  components.  Results  indicate  that  the  best 
localization  estimates  are  achieved  when  the  platform  is  on 
cardinal  headings,  and  when  the  -argst  moment  has  a  strong 
vertical   component. 
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I.    INTRODOCTIDH 

A.       BACKGROUND 

Ever  since  the  operational  introduction  of  the  MAD 
(Magnetic  Anomaly  Detection)  system  in  1944,  it  has  been  a 
single  sensor  system  limited  to  indicating  detection  only. 
The  short  ranges  of  the  aarly  versions  allowed  a  general 
localization  if  a  signal  was  detected.  The  use  of  multiple 
sensors  should  allow  for  a  good  estimate  of  the  target's 
location  relative  to  the  sensor  platform,  assuming  that  the 
arrangement  of  the  sensors  is  such  that  the  received  signals 
are  linearly  independent  of  one  another.  Previous  attempts 
to  solve  this  problem  have  in  general  focused  on  the  use  of 
only  two  sensors,  but  did  not  develop  a  usable  solution  to 
operational  completion,  ifork  has  been  dene  on  this  topic  by 
wynn  et  al. ,  [ Hef .  1],  using  a  superconducting  aradiometer 
array.  The  measured  signals  in  their  work  were  subjected  to 
a  novel  signal  processing  technique  which,  in  the  labora- 
tory, was  used  to  invert  the  dipole  field  equation  to  deter- 
mine both  the  position  and  the  moaent  vectors  of  a  dipole 
signal  source.  This  thesis  approaches  the  same  proolem  by 
simulating  five  standard  total  field  measurements.  As  shown 
in  Figure  1.1  ,  four  magnetometers  moving  in  parallel  paths 
are  used  to  obtain  these  measurements.  The  advantage  of  this 
approach  compared  to  that  of  Wyr.n  et  al.  is  that  it  utilizes 
existing  technology  in  conjunction  with  approximations, 
delineated  in  Chapter  3,  in  orier  to  produce  a  usable 
result. 


B.       APPROACH 


a 


Th3  problem  of  localizing  a  target  from  any  platfor 
involves  at  lsast  three  unknowns:  r„  ,  ry  ,  and  rt  ,  the  compo- 
nents of  the  position  vector  from  the  platform  to  the 
target.  In  the  case  of  using  HAD,  three  additional  unknowns 
mx,  m  ,  and  mT ,  the  components  of  the  target's  magnetic 
moment,  are  also  involved.  The  result  is  that  fchere  are  six 
unknowns  that  must  be  solved  for  in  order  to  localize  a 
target  using  MAD.  Since  a  normal  sized  magnetic  targe-1-  car- 
be  approximated  by  the  magnetic  field  of  a  dipole  [Ref.  2], 
it  was  decided  to  use  the  dipole  aquation  as  described  in 
[Ref.  3]  to  generate  the  signals  used  in  this  simulation. 
Equation  1.1  describes  the  MAD  signal,  and  will  be  developed 
from    the    dipoLe    equation    in    chapter   2.      Equation    1.1    is   r.on- 

S(r)=-(m^  +  me+m^  /r^+3 (a^* mr+rn^  (^-y^/rs         (eqn    1.1) 

linear  in  the  unknowns.  The  position  vector  in  equation  1. 1 
is  defined  as  going  from  the  target  to  the  platform;  this  is 
a  consequence  of  using  the  dipole  equation,  and  will  be 
explained  later. 

The  cartesian  coordinates  are  lefined  in  Figure  1.2  , 
The  z  axis  is  in  the  direction  of  platform  motion,  and  the  y 
axis  is  vertical  down.  The  components  of  the  earth's  field 
unit  vector,  ex  ,  ey ,  and  e2  ,  are  assumed  to  be  known  quanti- 
ties. Since  it  is  very  difficult  to  explicitly  solve  a  set 
of  non-linear  equations  in  six  unknowns,  the  problem  will  be 
simplified  by  making  assumptions  and  dividing  if  into  two 
special  cases.  The  first  case  is  a  trivial  ona  in  which  the 
target's  position  is  assumed  known.  For  this  case,  equation 
1.1  is  linear  in  the  three  unknowns  nx  ,  mr  ,  and  m?  .  If  the 
three  measurements  are  linearly  independent,  then  Cramer's 
Rule,    [Ref.    4],     may    be  used    to    solve    for    these   unknowns.    The 
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second  case  is  the  mere  difficult  dh=,  and  it  is  the  problem 
investigated  in  this  thesis.  In  this  case,  the  assumption  is 
made  that  the  targets  magnetic  moment  is  known.  That 
reduces  eguation  1.1  into  a  non-linear  equation  in  only 
three  unknowns,  rK  ,  ry,  and  r2 .  If  this  simplified  version 
of  equation  1.  1  could  be  manipulated  into  a  linear  form  by 
using  several  sensors  and  appropriate  combinations  of 
signals,  the  position  of  the  target  could  easily  be  solved 
for. 

This  thesis  will  develop  an  approximate  solution  to  -his 
problem  by  using  combinations  of  the  total  field  signals  of 
an    idealized  multiple  MAD   system. 
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Figure   1.1        Spatial   Relationship   of   Sensor  Positions. 
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II.    THEORY 

The  total  field  magnetometer  measures  the  magnitude  of 
the  vector  sum  of  the  ambient  magnetic  field  of  the  earth 
and  xhe  magnetic  field  of  a  magnetic  dipole  moment.  The 
field  of  the  dipole  will  used  to  represent  *rhe  magnetic 
field  of  the  target.  The  magnitude  of  the  vector  sum  is 
called   the   total   field  as   given   by 

rt  =  i  re  +  bA  i 

B*e  represents  the  earth's  magnetic  field  vector  and  b^ 
represents  the  field  due  to  a  magnetic  dipols.  In  the 
general  case  of  interest  here,  the  magnitude  of  the  dipole 
field  is  much  less  than  the  magnitude  of  the  sarth's  field. 
Under  certain  conditions,  a  useful  approximation  to  the 
total  field  may  be  derived  from  the  geometry  shown  in  Figure 
2.  1     . 


Figure    2.1        Vector   Relationship   of    B    ,    b^ ,    and   \. 

In  this  figure  a  very  small  vector  is  added  to  a  very  large 
one.  The  vectors  of  Figure  2.1  ars  not  to  scale.  b±  is  to 
be  considered  several  orders  of  magnitude  smaller  than  if,. 
The  applicatiDn  of  basic  geometric  formulas  for  the  sides  of 
a  triangle   and    the    angle    oetween    sides    yields 
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for  the  geometry  of  Figure  2.1  .  In  terms  cf  0  this  equation 
becomes 

Bt2    =    Be2    >   b^2   ♦    2bolBecose 

The  3e2  term  can  be  factored  out  of  the  right  hand  side  of 
the    eguation  to    produce 

Bt2    =    Be2{l    +     (bd/Be)2    f     2(b<JL/Be)cos9} 

Talcing   the   square  root  yields    equation   2.1    . 

Bt=Be£[U(b<A/Be)2*2(bjL/Bc)cos9]Vfe}  (eqn    2.1) 

Since  by  assumption  b^  is  much  lass  than  Be  ,  the  term 
(b^/Be)  2  is  much  less  than  2  (b^/Bg  )  cos9  and  can  be  consid- 
ered  negligible.    This  reduces   equation    2.1    to 

Bt»Bf  £[1  +  2(bJL/Be)  cos0]V*}  (eqn    2.2) 

Ths    general   expansion   formula   given    by 

(1+x)n     =    1  ♦    nx  ♦    n(n-1)x2/2!    ♦    Higher   Order  Terms 

can  be  applied  to  equation  2.1  .  For  the  case  here,  n=1/2, 
and  x=  (b^/Bg  )2 +  2  (b^/Bc)  cos9  .  Substituting  these  values  into 
the    general   expansion   formula    gives    aquation    2. 1    as 

Bt     =    Be{1    +    .5[  (b^/Be)  z    ♦    2  (bA/Be  )  cos9  ] 

~    1/8[  (ba/3e  )  2    +    (2b,L/3e)  cos9  ]2    ♦    H. O.T.J 

Squaring  the  appropriate  tsrm  and  multiplying  through  by  the 
coefficients  results    in 

3t    =   3e  {1    ♦    -5[(b-t/BB)2]    ♦    (b^/Beioos9   -    1/8[hgi/Be34 
-    U/S[  bA/Be  ]3cos9    -    4/8[  b^/Bg  ]2cos2  9    +    H.O.T.} 

The  cubed  and  higher  order  terms  will  be  considered  to  be 
negligible.  Dropping  those  terns  and  rearranging  the 
remaining   ones    yields  the   form 
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Bt    =   Be{1    «■    (b^/B^Jcose   -    .5(bA/Be)  2Cos29    ♦    .5(bol/=e)2} 

Using   the  relation    sin29   ♦    cos29   =    1,    this   reduces   tc 

Bt    =  Be  {1    ♦    (bA/Be)cos9  ♦    .5  (b^/Bg)  2sin29} 

After   multiplying     through    by      Be  ,      the      second   term      can   be 

— *     —* 
written     in  tar  as     of     the    dot      proiuct      ^ji*Be/3e  .      This     is 

reduced   to      bjL*e   if    the   unit      vector    e*  is    defined      as    Bg/Bg, 

and    the   resulting  eguation    is 

3t    =   Be    +    b^»e   ♦    .  5  (b^/Bft  )  2sin29 

Since  b±  is  on  the  order  of  1nT  while  3e  is  on  the  order  of 
50,000  nT,  the  error  term  involving  the  sine  is  negligible. 
For  the  stated  assumptions,  equation  2.1  can  be  represented 
by 

Bt    =   Be    +    b^»e  (eqn    2.3) 

In  this  thesis,  3e  will  be  assumed  constant.  In  addition, 
this  simulation  will  assume  that  no  noise  exists.  Thus,  the 
total  field  ai  a  g  re  to  meter  will  measure  the  anomaly  in  the 
earth's    field  defined   by 

S  (r)    =   b^  .2 

In  this  expression,  b^  represents  the  magnetic  field  of  a 
dipole.    b^    is   defined  [Ref.    3]   by   the    equation 

b^  (r)    =  -m/r3    +   3  (nUf)  r/r5  (eqn    2. a) 

m  is    the   magnetic  dipole    moment  in  units   of   nTft3.      r    is   the 

position   vector    to    the  field   point  and   is    measured    in    feet. 

If  b^  is  much  less  than  3€  ,  then  aquation  2.3  and  equation 
2.4    result   in  the  equation 
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S  (r)=-  (m»e)/r3+3  (m«r)  (r«e)/rs 


(aqn    2.5) 


This  equation  is  an  alternative  basis  for  the  derivation  of 
the  Anderson's  Functions,  found  in  [Ref.  5].  The  Anderson 
Functions  are  a  commonly  used  approximation  for  the  signal 
received   by   a  MAD   sensor.  The   computer   simulation    in   this 

thesis  is  based  on  equation  2.5  directly,  since  it  is  in  a 
representation  which  is  easier  tc  use  than  the  Anderson's 
function    representation   of    a   MAD    signal. 
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III.    NATORE    OP    THE    PSOBLEM 

This  thesis  will  develop  and  dam on strata  an  approxima- 
tion that  generates  a  reasonable  solution  of  the  target's 
position  relative  to  a  platform  bearing  several  "sensors"  to 
"measure"  the  magnetic  signal  as  the  platform  makes  a 
straight  line  encounter  near  the  target.  The  assumptions  are 
as    follows: 

1.)      The   dipola      equation   is   a   valid     representation  of 
the   magnetic    field    of   the    target. 
2.)    djl  is  several   orders    of   magnitude    smaller    than    3e . 
3.)       Measurements    are    done   using   only    straight    line    en- 
counters. 

4.)      The     target's   magnetic      dipole   moment      is    a      known 
vector. 
5.)    The   earth   magnetic    field   vector   is    known. 
6.)    The    gradient    of  the    earth's   field   is   zero. 
7.)    No   magnetic    ncise    of    any   kind     (environment,    sensor, 
etc)     exists. 

8.)         The  measurements      are    oriented      along   four      lines 
which    are     oriented  with    respect      to    the    flight      path   as 
shown    in    Figure    1.1. 
Since    the    dipole    equation   is    written      with    the   moment    at    the 
origin      of    a     coordinate   system,        the      calculations   of      the 
signal    field  values    used    in    the   simulation    are    made   with   the 
target    located    at   the  origin.      As    will    be   seen,       the    program 
performs   a    coordinate  transformation    in    order   to    produce  the 
output  of  target   position   relative      to   the   platform.         After 
this    coordinate    transformation,    the    origin    is   located   at   the 
sensor    platform    and    moves   with    it. 

The  five  measurements  are  made  along  four  separate 
flight    paths,    so   only   four    different    "sensors"    are    involved. 
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Sensor  4  could  be  replaced  by  inducing  a  time  delay  in  the 
measurement  of  sensor  1.  The  signals  generated  at  these 
measurement  points  are  measured  at  iach  time  step.  The  posi- 
tion vectors  can  all  be  described  relative  to  the  origin, 
located  at  the  target,  in  terms  of  the  position  vector  cf 
sensor  1.  The  measured  signals  are  represented  by  the 
following   equations: 

s  (f1)=-  { tmj&nff  wgif  (r1)  2-3  (m^c+i^y+m^i  (  xyyyz^J )  /r  1  s 
S<r2)*-{<Mr+Wl%)  (r2)  2-3  (ayc+n^y-dj+ayz)  (x^[y-d]yz^)  }/r25 
S  (r5)  =- {  (m^+rnye/ mep  (r3)  2-3  (  mj[x-d  ]  ♦  a^y  ♦  m^z )  ([  x-d  ]e+ye+  zej  }  /r 3s 
S(rt)=-  {(m^nyy-n^)  (rU)  *-3  ( nyc«-nyym^  z-d  ])  (xe/y  a+[  z-d  ]ej  }  /r  us 
S  (r5)=-  ((m^-ayy- mep  (r5)  2-3  (  mf  x+d  l+iy  +  m^z)  ([  x+d  ]e  +  ye+ze^  }  /r  55 

For  convenience,  the  following  dummy  variables  are  used  in 
the    computer  program    and   are   dsfinad    as: 

G    =    mKex  +  myey  ♦  m,^ 

V1  =    3mxe,(    -    G 

V2  =   3(mxey    +    myex ) 

V3  =    3myey    -    G 

74  =    3(mye2    ♦    m2ey  ) 

V5  =   3mze^   -    G 

76  =    3(mTex    -    mxeT  ) 

51  =    S  (r*l) 

52  =    S  (f"2) 

53  =    5(r3) 

54  =    S  (ft) 
35  =    S(r$) 

The  use  of  these  variables  reduces  the  above  equations  to 
the  following  form,  which  are  the  iifining  equations  for  the 
signals    in   the    simulation. 

S1=  (V1x2  +  V2xy+V3y2+  V4yz  +  V5  z2  + V6xz>  /r  1s 
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S2=  (71x2+7  2xy-7  2xd+73y2-27  3yd+73d2+7  4yz-74dz+75z2+76xz)/r25 
S3=  (V1x2-271xd+71d2+72xy-7  2dy+73y2*7  4yz+75z2+76xz-76dz)/r35 
S4  =  (71x2+7  2xy+7  3y2  +  74yz-74yd+75z2-2  75zd+75d2+76xz-76xd)/r4S 

S5=  (71x2+271xd+71d2+72xy+7  2dy+73y2+7  4yz+7  5z2+76xz+76dz)/r55 

The  unknowns  of  interest  are  x,  y,  and  z,  which  are  related 
to    the   position    vectors   as    follows: 

r"1  ■  xt  +    yf  ♦   z1c  (eqn  3.1) 

r5  =  xt   +    (y-d)t  +   zk  (eqn  3.2) 

r5  =  (x-d)t  +    ff  +    zfc  (eqn  3.3) 

rt*  -  x"?  +    y^  ♦    (z-d)lc  (eqn  3.4) 

r$  =  (x+d)"5  +    yj*  ♦    zlc  (eqn  3.5) 

To  linearize  the  problem  it  will  be  assumed  tha*  the 
denominators  of  the  equations  for  SI,  S2,  53,  34,  and  S5 
given   above   are    all    equal   to   r,    i.e., 

r    =   r1    =   r2  =    r3    =    r4    =   r5 
Tha    equations   for   S1r    52,    S3,    54,      and    S5    with  the    denomina- 
tors   equal   can    then    be  subtracted      froa    one   another    to    yield 
the    following: 

S2-S3=  {(72-27  1)  x*  (273-7  2)  y+ (74-76)  z+ ( 7 1 -73)  1}  (-d/rS) 
S2-S5={(72+271)  x  +  (273+72)y  +  (74+75)  z  +  (71-73)  i}  (-i/r5) 
S3-S1=  {271x  +  72y  +  76z-71d}  (-d/r^) 
S5-S1={-271x-72y-76z-71d}  (-d/r5» 

S  3-5  4=  {(271 -7 6)  x  +  (72-74)  y+  (76-275  )  z  +  (75-71  )  d}  (-d/r^) 
S5-S4=  {-  (27  1+76)  x- (72  +  7  4  )  y-  (76+27  5)  z+  (75-71)  1}   (-d/rS) 

Since  the  denominators  are  assumed  to  be  identical,  if  these 
six  equations  are  divided  in  pairs  the  coefficients 
involving  the  non-linear  r5  term  and  the  quantity  d  are 
canceled  out.  This  results  in  tha  following  equations  -hat 
define   the    new    variables    A1 ,    A2,    and    A3    as 

A1     =     (S2-S3)  /  (S2-S5)  (eqn    3.5) 

A2    =     (S3-S1)  /  (S5-S1)  (eqn    3.7) 

A3    =    (S3-54)  /(S5-S4)  (eqn    3.3) 
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Defining   the  following  variables    for    convenience, 

C1  =  A1  (72+271)    -     (V2-2V1  ) 

C2  =  A1  (273+72)    -     (2V3-72) 

C3  -  A1  (74+76)     -     (74-76) 

K1  =  [  (71-73)     -    A1  (V1-V3)    ]d 

C<4  =  A2  (2V1)     +    2V1 

C5  =  *2(72)    +    V2 

C6  =  \2  (V6)    +    V6 

K2  =  [  V1    -    A2  (V1)  ]d 

C7  =  A3  (76+271)    -     (V6-2V1) 

C8  =  A3  (74+72)     -     (VU-72) 

C9  =  A3  (275+76)    -     (275-76) 

K3  =  [  (71-75)     -    A3  (71-75)    ]d 

equations   3.6   -    3.3    can    be    writ-en   in    the    following    form   as 


C1x    ♦   C2y    +   C3z   =    K1 

C4x    +    C5y    +  C6z   =    K2 
C7x    +    C8y    +    C9z   =    K3 


(eqn  3.9) 
(eqn  3.  10) 
(eqn    3.  11) 

These  three  equations  are  linear  in  terms  of  x,  y,  and  z.  In 
addition,  they  are  linearly  independent  as  shown  by  using 
the  program  in  Appendix  3.  A  sample  output  demonstrating  the 
linear  independence  is  shown  in  Appendix  C.  Since  they  are 
linearly  independent,  Cramer's  RuLa  can  be  used  to  solve 
equations  3.9  -  3.11  for  x,  y,  ani  z.  The  computer  program 
in  Appendix  A  was  written  in  terms  of  tnese  aquations.  The 
program  assumes  that  the  magnitude  of  the  target's  magnetic 
moment  is  5x198  nTft3.  A  representative  moment  was  selected 
by  using  the  values  listed  by  Fromm  in  [Ref.  2],  who  lists  a 
moment  of  1 08  to  2x10a  cgs  units  fsr  a  submarine.  One  cgs 
unit  is  approximately  3.35  nTft3.  The  value  used  in  the 
program  is  3.  rounded  out  average  of  Fromm's  values. 
Computations  are  done  using  the  standard  convention  with  the 
target   dipole  at   the   origin.      However,    as    indicated   earlier, 
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the  output  is  corrected  by  a  coordinate  transformation  so 
that  the  platform  is  at  the  origin,  and  the  target  position 
is   given    with  respect  to    the   platform. 


IV.     DISCUSSION 

A.       GENERAL 

The  program  in  Appendix  B  was  written  to  check  whether 
the  results  generated  by  the  approximation  aquations  were 
reasonably  close  to  the  actual  coordinates  of  the  target. 
The  program  was  modified  into  the  form  shewn  in  Appendix  A 
to  graphically  provide  a  comparison  between  the  calculated 
and  the  actual  position  coordinates.  Initial  parameters  for 
the  knowns  in  the  equations  were  a  platform  heading  of  30° 
magnetic,  an  earth  field  vector  of  70 o  down  from  the  hori- 
zontal, an d  a  target  dipole  moment  vector  with  a  vertical 
component  50°  down  from  the  horizontal  and  a  horizontal 
component  oriented  355°  from  magnetic  north.  These  initial 
parameters  have  no  special  significance,  and  merely  repre- 
sent a  convenient  starting  point.  As  discussed  earlier,  the 
magnitude  of  the  target  dipoie  moment  was  set  at  the  value 
of  5X108  r.Tft3.  The  spacing  between  measurement  positions, 
d,  was  taken  as  50  feet,  and  the  platform  passed  directly 
overhead  the  target  at  1000  faet.  This  and  ail  subsequent 
runs  were  started  at  5000  feet  prior  to  CPa  (closes*  point 
of  approach)  with  the  platform  advancing  50  feet  per  time 
step.  One  position  calculation  was  completed  at  each  time 
step,  the  simulation  ceasing  when  the  platform  was  5000  feet 
past  CPA.  For  interpretation  purposes,  calculated  target 
positions  are  shown  as  the  sensor  platform  moves  10,000  feet 
past  the  target  in  a  straight  line  encounter,  with  C?A 
occurring  at  the  5000  foot  point.  Neither  the  CPA  nor  the 
target's  location  are  in  any  way  considered  known  in  the 
localization  process.  However,  as  the  area  of  interest  is 
the    behavior      of   the    simulation      close    to    CPA,         the   closest 
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Pigure  4.1        Initial   Run,    x   position, 


point      of      approach      is   always      centarsd      in     ths      following 
figures. 
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Figure  4.2   Initial  Run,  y  position. 

From  the  orientation  of  the  coordinate  system  shown  in 
Figure  1.2  ,  z  is  in  the  direction  of  flight,  y  is  straight 
down,   and  x  is  off  the  right  side  of   the  platform.   Since 
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Figure  4.3   Initial  Runr  z  position. 

this  thesis  investigates  straight  line  encounters  only,  the 
x(actual)  and  y (actual)  coordinates  will  remain  constant 
throughout  any  run  for  the  coordinate  orientation  used.   The 
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Figure  4.4        Initial    Run,   corrected   y  position, 

five  measurements  will  therefore  be  made  with  their  respec- 
tive x  and  y  coordinates  as  constants.  The  z  coordinate  will 
vary.        This   means      that      for      the    initial      pass      inscribed. 
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x (actual)  would  be  constant  at  z  =  ro,  y  (actual)  would  be 
constant  at  1300  feet,  while  z  (actual)  would  progress  from 
5000  feet  to  -5000  feet.  Figures  4.1  -  4.3  show  a  compar- 
ison of  the  actual  and  calculated  target  positions  as  the 
platform  compLetes  the  initial  run  described  earlier.  It 
will  be  noted  that  the  calculated  and  actual  x,  y,  and  z 
positions  are  almost  identical  at  the  start  of  the  run.  This 
is  a  coincidence.  "Crossover"  points  such  as  this  occur  for 
all  three  coordinates  simultaneously  during  a  run,  failing 
closer  or  further  away  from  C?A,  depending  upon  the  choice 
of   parameters. 

In  this  initial  run,  it  is  obvious  that  the  "worst" 
estimates  occur  for  the  y  coordinate.  This  is  less  signifi- 
cant than  it  seems  at  first  glancs,  however,  as  the  y  coor- 
dinate is  the  vertical  separation  between  the  platform  and 
the  target.  If  the  target  is  submerged,  this  value  is 
unknown,  but  the  altitude  of  the  platform  is  known.  The 
target  can  be  no  closer  to  the  platform  than  this  altitude, 
since  the  target  is  physically  unable  to  rise  above  the 
ocean  surface.  As  an  example,  if  a  target  depth  of  200  feat 
is  assumed,  the  initial  y  separation  of  1000  feet  forces  the 
altitude  in  this  simulation  to  be  900  feet.  Using  this  value 
as  a  bound  on  y ,  it  can  be  seen  in  Figure  4.4  that  the  y 
calculated  position  is  somewhat  improved.  Accordingly,  it 
seems  that  the  particular  arrangements  of  sensors  shown  in 
Figure  1.1  does  show  promise  in  the  problem  of  utilizing 
this  geometry  of  Multiple  MAD  for  accurate  localization 
purposes.  Accordingly,  variations  of  parameters  were  tried 
in  order  to  determine  what  effect  changing  them  would  have. 
In  all  cases,  parameters  not  specifically  mentioned  in  a 
section  will  retain  the  values  as  for  the  initial  run 
described  above.  The  singularity-like  behavior  occur ing  at 
various  places  on  the  graphs  are  a  result  of  the  determinant 
of  the  denominator  in  the  Cramer's  Rule  subroutine  having  a 
value   close    to    zero. 
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B.        HEADING    CHANGES 

In  this  section,  the  localization  process  will  be  inves- 
tigated as  a  function  of  heading.  Compared  to  the  initial 
conditions  described  above,  the  crossover  points  were  found 
to  be  moved  further  our  past  5000  feet  prior  to  CPA  on  a 
heading  of  0°.  While  a  specific  example  of  this  is  not 
shown,  the  behavior  of  the  crossover  points  for  the  platform 
headings  of  30°  and  60°  may  be  observed  by  comparing  Figures 
4.1  -  4.4  to  the  respective  portions  of  Figure  4.5  .  The  x 
position  estimates  were  better  for  a  North-South  heading 
than  for  any  other  platform  heading,  although  the  estimates 
on  East-West  headings  were  also  very  good.  For  a  90°  heading 
there  are  two  crossover  points  located  equidistant  at  1200 
feet  from  CPA.  Calculated  positions  en  all  headings  are 
generally  gcod,  excepting  those  around  60°,  where  the  x 
estimates  shows  a  large  discrepancy  at  distances  greater 
than  3000  feet  prior  to  CPA.  The  y  5stimar.es  are  not  signif- 
icantly affected,  but  the  z  estimates  are  good  only  within  ± 
1200  feet  of  CPA.  Nevertheless,  even  in  this  worst-case 
situation,  reasonable  calculated  coordinates  are  available 
for    localization   along  a    significant    portion    of   the   track. 

As  shown  in  Figures  4.5  -  4.8,  headings  in  the  quadrant 
0o  -  900  reasonably  describe  what  happens  on  any  heading,  as 
the  graphs  for  analogous  headinas  (50o,  120°,  24QQ,  and  300° 
in  this  example)  have  similar  shapes,  although  the  orienta- 
tion may  differ  between  quadrants.  These  similarities  are 
apparently  due  to  the  fact  that  analogous  heading  in  quad- 
rants 1  and  3  and  quadrants  2  and  4  are  in  fact  reciprocal 
headings,  and  so  one  would  expect  the  graphs  to  maintain  the 
same  shape;  however,  the  ending  point  of  one  would  represent 
the  starting  point  of  the  other.  The  relationship  between 
adjacent   quadrants    probably    exists   for   similar  reasons. 
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C.  TARGET    DIPOLE   MOMENT    CHANGES 

In  this  section,  the  horizontal  component  of  m  was 
oriented  in  the  0°  (magnetic  North)  direction  vice  355°,  and 
the  vertical  component  of  m  was  varied.  The  best  positioning 
was  obtained  for  m  oriented  vertically  up, i.e.,  no  hori- 
zontal component,  and  an  example  for  that  case  is  shown  in 
figure  4.9  .  Poor  positioning  was  obtained  for  the  vertical 
component  of  a  pointing  between  the  horizontal  and  45°  up. 
The  worst  case  shown  is  Figure  4.10,  where  a  is  oriented 
horizontally  pointing  due  north,  with  the  moment  oriented  at 
an  angle  below  the  horizontal,  reasonable  positioning  esti- 
mates were  obtained.  There  was  Less  variaoility  in  posi- 
tioning for  moments  pointing  below  the  horizontal  than  for 
those   pointing    above   the   horizontal. 

The  results  for  variations  in  the  orientation  of  the 
target  dipole  moment  suggests  that  the  positioning  estimates 
are  less  sensitive  to  moment  orientation  changes  in  the 
horizontal  plane  as  compared  to  taa  vertical  plane.  For 
example,  by  keeping  the  vertical  component  of  1*  constant  at 
an  angle  of  53°  down  while  varying  the  horizontal  component 
direction  it  was  found  that  the  horizontal  component 
affected  the  positioning  much  less  than  variations  in  the 
vertical.  Illustrated  in  Figure  4.11,  the  x  estimates  of 
position  improve  as  the  horizontal  component  approaches  the 
platform  heading,  while  the  y  estimates  improve  as  the  hori- 
zontal component  approaches  alignment  with  magnetic  North. 
The    latter   case    is    not  shown. 

D.  CHANGES    IN    SENSOB    SPACING 

Changing  the  inter-sensor  spacing  had  at  most  a  weak 
effect  upon  the  position  output,  even  with  spacing  as  short 
as  5  feet  or  as  long  as  200  feet.  This  would  be  expected  in 
a  computer      simulation  where   sufficient    precision      is 
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available.  For  a  spacing  as  graat  as  1DD0  feet  sone  effect 
became  notioable.  in  the  real  world,  sensor  spacing  w2uld 
influence  *:he    localization   estimates. 

E.  ALTITUDE  CHANGES 

For  a  straight  overhead  pass,  positioning  information 
appaars  best  if  a  pass  with  a  500  to  1000  foot  separation 
between  the  target  and  the  platform  is  made  at  CPA.  At 
greater  altitudes,  the  accuracy  suffers  slightly  with  the 
increase  in  altitude,  while  at  lesser  ranges,  usable  posi- 
tioning information  is  obtained  only  relatively  close  to 
CPA.  For  example,  for  a  vertical  separation  of  400  feet, 
positioning  data  was  very  inaccurate  at  slant  ranges  greater 
than    100  0   feet. 

F.  NON-OVERHEAD    PASSES 

This  program  gives  poor  results-  for  passes  that  are  no* 
directly  over  or  directly  off  to  one  side  of  the  target.  The 
worst  situation  observed  in  any  variation  is  illustrated  in 
Figure  4.12  .  To  generate  this  figure,  a  pass  was  made  such 
that  the  actual  x  coordinate  was  approximately  equal  to  the 
actual  y  coordinate.  The  best  positioning  was  achieved  by 
overhead  passes,  but  low  off  to  the  side  passes  can  yield 
good    results  as    well. 

6.       AVERAGES 

The  averages  resulting  from  a  number  of  runs  are 
presented   in  Table    I    .  These   selected   averages,      covering 

all  of  the  variations  discussed  in  this  chapter,  show  little 
consistency.  This  is  probably  due  in  part  to  the  existence 
of  the  singularities.  The  variability  is  also  a  function  of 
the    limits   of   travel    for      the   runs.       Nevertheless,      for   most 
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Figure    4.5        Platform   Heading   60°    Magnetic. 
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Figure    4.7        Platform   Heading   240<>    Magnetic. 
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Figurs    4.8        Platform    Heading    300°    Magnetic. 
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Figure  4.10   m  Horizontal  and  Orisntad  to  Magnetic  North. 
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cases,  -hey  do  show  reasonable  values,  and  by  themselves  io 
provide  limited  localization  data.  The  problem  fcr  inflight 
analysis  would  be  to  determine  the  bounds  of  a  run  in  -he 
absence  of   any    other   target    positioning    data. 

H.       IHPOBTANCE    OF   C0NSIST2NCY 

An  important  observation  resulting  from  the  analysis  of 
the  graphs,  including  those  not  shown  as  examples  in  this 
chapter,  is  that  the  rate  of  change  of  the  calculated  coor- 
dinates gives  some  indication  of  the  validity  of  -he  posi- 
tioning data.  If  the  x  and  y  graphs  have  a  slope  close  to 
zero,  and  if  the  z  graph  has  a  slope  of  -45°  for  the  calcu- 
lated coordinate,  then  positioning  data  is  likely  to  be 
comparatively  accurate.  If  all  three  graphs  fulfill  this 
condition  simultaneously,  then  positioning  data  appears  to 
be    usable. 

I.        IMPLEMENTATION 

The  possible  implementation  of  the  process  described  in 
tnis  thesis  *ould  be  to  assume  a  magnetic  moment  for  the 
target,  then  run  the  signal  data  through  the  processing  to 
generate  positioning  data,  and  then  compare  zhas^  results  to 
those  obtained  by  using  the  assumed  moment  to  generate  a 
signal  via  this  simulation.  If  the  agreement  between  results 
were  good,  one  could  assume  that  one  had  assumed  the  correct 
target  moment;  if  not,  one  would  assume  a  different  moment, 
and  gc  through  the  process  again,  until  one  had  a  reasonable 
target  moment,  which  would  be  vaiii  at  least  until  the 
target  changed  course.  (Jsing  a  computer  these  comparisons 
could  be  made  reasonably  quickly,  so  this  process  might  have 
operational    possibilities. 
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▼•    COHCLOSIOH 

For  the  assumptions  made,  the  positioning  estimates  are 
reasonable.  Also,  the  averages  taken  of  the  three  coordi- 
nates as  calculated  throughout  the  various  runs  show  that 
tha  x  calculation  averages  are  ia  general  at  worst  within 
500  feet  of  the  actual  x  value,  while  the  worst  y  calcula- 
tion averages  tend  to  be  within  1000  feet  of  the  actual  y 
value.  The  best  calculated  x  and  y  values  were  within  100 
and    200    feet,   respectively,    of  the   true   values. 

Trends   observed    were   as    follows: 

1.  It  is  best  to  fly  the  sensor  platform  directly  over 
the  target,  or  alternatively  to  have  a  'distant' 
(between  1000  and  2000  foot)  CPA  off  to  one  side  while 
minimizing  the  target  to  platform  separation  in  the  y 
direction.  The  worst  situation  is  when  the  target  is  at 
a  point  that,  at  CPA,  has  a  depression  angle  of  45°  with 
respect   to  the   platform. 

2.  Best  vertical  separations  are  those  that  will  yield  a 
500  tc  1000  foot  vertical  distance  between  target  an  1 
platform.  Short  vertical  separations  should  be  avoided 
if    good   localization  is    desired. 

3.  Bes*  results  are  obtained  if  headings  are  northerly 
or  southerly.  Poor  results  are  achieved  wi*-.h  headings 
close  to  600,  120O  ,  240°,  and  303°.  East-tfest  headings, 
while  not  as  good  as  North-South,  do  provide  good  posi- 
tioning  data    close    to   CPA. 

4.  Good  results  may  be  expected  if  there  is  a  large 
vertical  moment,  especially  if  it  is  vertically  up. 
Poor  positioning  results  are  obtained  with  this  simula- 
tion  if  the   moment    has    a    large    horizontal   component. 
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5.      Results  are   slightly    better      if    the    target    moment   is 

lined   up      with   either      the   axis      of   flight      or    with      the 

earth's   field,      although    this      effect   is    often    dominated 

by    the   more  gross    effects    caused  by    the    ether   variables. 

A      suggested   approach      for   further      study      is    to      inject 

noise   into   the      simulation    in   order   to    see    how      robust    it    is 

in      the   accuracy      of      the      position   calculations.  If      the 

results   remain    reasonable,      this   approximation   or    an    altered 

form    of    it   may    one    day  have    real-world   applicability. 
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APPENDIX  A 


COMPUTER  PROGRAM  FOR  PLOT  OUTPUT 


oooooooooooooooooooooooooooooooooooooooooooooooo 
^ovjrn>ru^<)r^co(>o»^og<r»vrir\^or^aoo*0'^<N<^>*'if\^r^aocFO^cvjrn  or*- co 

oooooooooooooooooooooooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooooooooooooooooooooooo 


ooooo 

o 

o 

o 

o 

o 

o 

o 

O      Z2T 

o 
o 
o 
o 
o 
o 
o 
u 
o 
o 
o 
o 
o 
o 
o 
o 
u 
o 
o 
o 
o 
u 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
u 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
ooooo 


<QdOO 

i-.ua: 

1-OU.O 
LU)-o 

a: -jo 
o     < 

UJCL 

JIO 

<xt- 

o      • 

t—  00 1— 

o<lu 
woo; 

aoH 
I— I— 'JJ 

"3  — X 
CLoOt— 
l-O 
DQ.I- 
O      00 

_)<Q. 

<~ 

Ooooo 
— «LUOO 
IH< 

a.  a:  a. 
<<I 

QCOJ 
O       LU 

Q> 

-dLULU 
H-_ I 
00< 

oo< 
ou 
o<>- 
adox 
a.     o 

sz< 

«.aL 

OOOOO 

az 

0dO< 

t— oo 

00»— LU 
i— lOO^ 

XL_>< 

i-cls: 


2:2:2:2:2:2:2:2:2: 

uuu      » 

o>-  * 

l_)LUr\j 

ut-  « 

o  -*: 

oz  • 

OO0> 

o<o 

Ooi  -  — O 
O  »J0  r\j<_) 
0*-«0 

OQ.    » 

O   •J>" 

uoo 

o  -  * 

omo 

OoOO 
0<\|    •> 

Ooo^n 
o  -o 

Om  •■     — o 
<JoO>3"      s:r\i 
OMO 
OoO    * 
O  »<-<*> 

uoo 

Odd   - 
O   »<N 
O00O<"-'O 
O   *  *t— 0<N 

OQhujZ- 

o  ••ocOGdrNi 

UN    *Qd    "O 

o  ►<,o  *—  «. 
ox<m-» 

O  *  "U-  *«■  -^ 

OXf\i_JLOO 

u  •*<<   *<M 

on  -a:—.— 

•     ouj'-<  »^r\> 

5£      O   »<IvO  — o 

O       0>    *>N    •» 

<  olunT  •»  ~— 
ac     o  •»</>m-»""* 

►—    oxro>ino 

oujoo  •»— <\J 

o    o  ••  •%*•>•>-. 

Z       ON^>    »X 

a    usw  •— -o 

<l       O>o0>— — 

us;  »  •..>*  — 

h-  O  ».*^(\|  O 
LU  OXoO>2CvJ 
LU       O^CO    •O  — 

lu     u     oo— *>— nvi 

OoO    »>OOct 
O      O*  LU   -2   »0 


xsjtzssssxsssssrssrsxjrszsrxssrrrs: 


a:o 
o> 

Xrr> 

a:  — 

"CD 


o  » 
-•< 

LUO 


LO-^ 


Zee 


LO.^ 


o_)aj<*-yj— . 


O  0<XoO 
LT\  OLUQ.iT\»-'0 
OQCoCOOQCNi 
h-  O  **  ♦ 
<       O 

o 
o 

OOO 


LUO 

'JO 
aru_ 

<  I 
>-< 

OS 

2< 
<o 

OO. 
•-0 
o 

LULU 

>o 

o 

CI— 

-J1-* 

LU2  • 
— «OLU 

u.<al 

3IUJ 
O      X 

(— LUO 
LUOOUI 

J3<LU 

<  t— 
srz 

t— LU 
OO"— • 

-  3!> 

X      _l 

—  »UJ 

CCCC> 

<o«— 

LUh-t— 

oo 

o>ai 

2      LU 

<LU« 

LU  — 
XU.LU 

s:o< 
a:— 

oi-  - 

U.LUQ 
I— ZLU 
<003 

axo 


00 
UJ 
LU 

ac 
o 

LU 

— 

o 
o 

I 
o 


O 

LL 

t— 
< 

U 

— 
O 

o 


< 

LU 

X 


o 
o 


o 

—I 

LU 

■a. 
o 


'XI 

o 
< 

3! 


< 


LL 

o 

•JJ 

— 
o 

cz 
<I 

a. 

** 

a 


o 
a 

ii 

+ 


2 

LU 

a 

2: 


LU 

O 
< 

ST 


LU 

o 


a 

- 
o 

< 

a. 

o 


2 

LU 

ST 

a 
o 

LU 

O 
< 


LU 
O 

< 


U- 

o 


o- 


000 

LU< 

oclu 
O  II 

-10 


o  » 

NX 

— it- 
ad  ex: 
oo 

iz 


cca     oca. 

LU<^       LU>— 1 


OCT. 
uua 


o 

a 


t— LU     t— o     t— »-» 


a£2Q  <<< 
3<  «*cx:a:a: 
X^-O 

<co<#  •«■  * 

1  o-<>- 

LULU<5< 


ad  11  2: 

LU       2C 

l-OO 

-»zi-*zai-ZQ-»z     j-a. 

ouj     -^uj     osiuj     colu  »     -•□2a#'siuju.i- 

lo     -*lo     -»m     -»Ln     h->->o<Ql<         ^-X_juj 

—    #     ••--»■     «••    *     ►•    OO*  LULU1— LUO'—       Q.<03 

vO—   »-v0*-   ».sO—    .^O—lu   •>      XOXOQl  II 

*-.       LO— .       LOw       LO^-        3LO  II  t  |t    ||    || 

H-»-.      I— ^»      1—^-      I—       *-        II    II    II  -t  II  < 

LU<       LU<       LU<       LU<  II         <  I— LU<< 

t-Snt-SOt- 5-Qt—  S:       Ot-LU<<t  ||  2LULUU.I— 

— cx:<'--icx:<«-'cr.<t>-icx:o<LULULL»—     oxx-juj 

JdO  LUQdOLUCX:  O  UJQC  O  O^  LUX  X  — JLU>-«  <  t—  a.  <  ca 

3u_cx:3u_a:2LLcr:3:u.  I  act—  q.<.cqq.ccoc<x.cccc 

o         -—        cm         m 
lt>         m        in         tr\ 


- 
00 

X 

t- 

ad 
< 

LU 

'JJ 


< 


LU2 

2:— 

00 

so. 

'JJOO 

— ,— 

ox 

Q-i- 

— 

LU 

IQ 


LU< 

O-l 

o 
000 
I— —I 
2< 

LUO 

X 

OUJ 
a-od 
21<l 

O 
Oct 

o 

Nl- 

*o 

>-LU 

-•> 

X 

H 

LUi-i 

J.2 

l-O 


LU 

CO 
ad 


oO 

o 


u_— 1 
<i< 

QdU. 
— -J 
oo< 

Oad 
o— 
02 
«— »i— 1 

■»■  00 

— o 
o— 

Q» 

«o 

LOO 
I       • 

—»m 

11  11 

x> 
2:2: 


ooooo 


0000 


44 


oooooooooooooooooooooooooooooooooooooooooooooooo 

^-  in  m  m  m  m  m  m  m  m  in  ^o  so  >o  o  «o  o  -o  vO  ^o  o  r-  Is-  [•»  Is-  f*-  r*-  r*-  r-  p-  p-  ao  co  ao  co  ao  co  ao  oo  00000^000^00* 

oooooooooooooooooooooooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooooooooooooocooooooooo 
ooooooooooocoooooooooooooooooooooooooooooooooooo 


O 

0 

XX 

LU 

— < 

> 

Utf 

■— • 

<to 

• 

LU 

a.a 

UJ 

O 

<S>CC 

cc 

UJ 

a. 

LU 

C£    • 

QC 

X 

h- 

ouj 

'JJZ 

oox 

Q 

50»-< 

x>- 

LU 

O 

UJ 

h- 

xa. 

OO    •» 

< 

»— 

a 

—I 

00 

LULU 

O 

•^-4 

Xa: 

O 

ZI 

-* 

>-•-• 

_l 

x>- 

— » 

to 

< 

1— 

«NJ 

GOLU 

O 

H-w 

* 

<Q 

LU 

<< 

X 

* 

_IGO 

CL 

ac 

—) 

-J—* 

< 

LU'U 

<m0 

-^ 

UJ 

t— 

rvl 

*"• 

35  LU 

m 

LU< 

>— 

< 

O 

-l_l 

h— ■» 

WZ 

1 

XX 

+ 

lu— . 

** 

<< 

0.O 

01  uu 

■»• 

X 

f\i 

—     1 

— . 

aruu 

LU 

•"O 

o< 

(\|^ 

*-x 

UJ 

r-t 

00 

0 

-»•  rsj 

oOO. 

X 

— » 

—» 

— « 

<r 

ct 

LU 

it-'XJ 

a  a: 

O. 

— * 

-^ 

— » 

a: 

O 

Xoo 

—  * 

o«- 

a:~» 

X 

> 

X 

XLL>- 

00 

h~—t 

-5  IS 

QZ 

>-M 

UJ 

UJ 

UJ 

00 1—2 

X 

—  3! 

**ii« 

00  UJ 

■»■ 

* 

# 

X      >- 

<u 

Zoo 

>w 

#00 

o# 

> 

rS 

M 

lu  1  a 

00 

xa: 

^-  + 

— Q 

OM 

X 

X 

X 

OO       Q. 

0 

^^^■^ 

QS 

»■» 

«- « 

— » 

UJ 

LL 

O00 

+  > 

<l* 

— O+O+O+ 

u_a£oo 

O 

LUZ 

LU 

u.-^ 

#  + 

—» 

— » 

« 

Olu« 

oOLU 

— - «■ 

_i-» 

—fm   | 

>•  1 

rvl  | 

M 

XX 

00 

<00 

CNJ> 

<i< 

-»>- 

UJ 

LU 

LU 

z     t- 

X 

ca 

*z 

Oil 

«.<LU  — 

•■*-. 

'♦  — 

* 

Oh- 

O 

LU 

-*  — 

«- LU<t>—  ■*  XX»NX 

«oik 

O      OO 

— 

•»:> 

—  + 

oox»-lu>-ujxluxlux 

1— oO< 

OOO 

0  »o  •  • 

H 

QO 

a 

0 

_j« 

-9— » 

OI-UJIS# 

*— * 

—  * 

m 

M_ | 

QQQ 

•  0  »oo 

►"^ 

<u_ 

—  X 

oorxi— — x 

— > 

—IS 

1- 

00  GOO 

•    •    • 

00000 

00 

1   1 

+ 

1 

X 

XUJ 

o— h-cc  +  s    2:    2: 

Q.-.UJ 

a.     ex 
s  uu 

OOO 

00 

000000 

•  OOOO"^ 

om— «m-*  1 

U 

u. 

OLU 

*-* 

•1— X 

#a- 

-cOXo 

1    0 

1    a 

1 

1  I.I  II  ll  1 

— 1,— __l 

—  — .— 

ll   ll   1 

OO  h- 

00; 

rOLUQOOOQO 
"O*     »Q    »Q    .Q 

-JQt- 

<r  _i 

-*mO 
r  1  0  11  11  11  11  11  11 

-J 
< 

-IU. 

K  0 

□  Q« 

X  >  X  >  X  >  X  >"xlMpxjivl 

twi/ii-jxm 

1    •en 

•  (T 

1    • 

1-1     < 

• 

1— • 

<o 

-^00  II 

m  1  a~x~ 

•m«- 

•rr»»- 

•CO 

H-LU 

h- 

II   II    II   II 

II  II  II 

II  II  II 

II  II 

X 

0 

«-• 

>-«OLU 

1 

n  — <  —»  cm  r\j  m  ro 

— 

h-X 

O       UJ 

11  11  1 

1   It        11 

11  it 

II  II 

II 

^»XCQ 

X 

00 

CM  II  2! 

II 

•"••si 

_t<M_i 

II  CQ<CQ<CQ< 

►— • 

<nj  t\j  corn  sj- 

<< 

t— 

rvlX^M      '-l-  > 

irxim^MTv  vTi 

—  *— 

OOOOOO 
3UUOUOU 

LU 

LU 

OOO 

qq:q 

Slululuo>;>;>:>;>> 

LUGO  GO 

X>ivli 

LU 

X>X>X>X>i>JPnIMI~nI 

J_—'-J 

— 

J->— 

h-Qi 

h- 

0 

ro 

H-<1 

00000 


OOO 


0000 


45 


oooooooooooooooooooooooooooooooooooooooooooooooo 

r»»coo»Ot-<r\jm^mv0^coo»0^c\jm>*-mNOr»»<DOO-*c\jmv^m^r^co 

O*  0*0*0000  0000  00 -*-*-*-*  »-««-<'-4-*— < -^rgcN<Nl(\i  eg  <\j  ro  f^fvjfsirn  corner  rorOrnroro  rnvj-vrNl-^"  >4" 

oooooooooooooooooooooooooooooooooooooooooooooooo 
000000000000000000000000000000000000000000000000 


UJUJ 

1— CD 

< 

zo 

—a 

O 

-JLU 

l— 

UJ-5! 

3 

Q£ 

Quo 

UJ 

1—00 

O 

< 

QC 

OC 

a  • 

UJOO 

oc 

o  — 

zo 

QC 

— i- 

act: 

< 

a 

zz 

Z^ 

_  — 

~< 

*s: 

*^-^ 

— .            M* 

z 

<o 

—cm 

—     o 

»•— < 

l-Z 

-<> 

—t     > 

zs 

UJ 

>  1 

>     1 

— i«^ 

ujO 

ujo 

♦  — 

•«■     — 

—  tvl 

iiZ 

QC 

ocn 

o     in 

->LU 

<LU 

<UJ 

o> 

Q     > 

—  * 

l-Q 

X 

♦*■ 

•     * 

rsl- 

OOI— 

(MO        O 

CM       OO 

■»■  -> 

LUUJ 

UJ 

—  Q 

—       O 

rvl«— 

OCX 

oz 

1     •      * 

1           •* 

SIM 

<i-  • 

z-» 

(N<N 

sO       CM 

w^. 

Z 

UJ 

>—      — 

— » 

^>      ^.-, 

+  + 

_i  »a 

cc  3: 

**M          — 

— 1 

■— '      »— — . 

— »— . 

^OO'— • 

LUQC 

__ 

> 

— .      — . 

—  :>- 

z^-*- 

U_UJ 

1    1    Ofl 

* 

1  cm  1  in 

-5UU 

Oac< 

U_l— 

:>:> 

O 

:>     > 

•-«->  —+                    —• 

«-»aiaE 

i— * 

— •— »  1     1 

o 

>— CMUJCM      — 

OOI— — 

act: 

i/i 

— .v*--. 

oO 

•       o 

oo 

—>*•—* 

*  -}■*■  ccs: — > 

X 

UJ 

i— 

—  C\J>> 

1— 

CM 

1— 

—  >vO> 

>■—  *H-C£ir\w 

oao 

UJO 

z 

— <;>>__ 

z 

<-.        * 

Z 

-^w>^* 

S>(JOLU#  S 

UJUJCX 

>C£ 

UJ 

>  +      * 

UJ 

UJ 

>     +* 

~wa:aH.#  o 

>QQ- 

oo 

— 

#  —  1  -* 

"-* 

+          — » 

t-4 

•»•  I  —co 

+  +*-*     e>-^ 

~cca. 

co  1 

o 

om    < 

o 

-» 

o 

O     in< 

—  #  Qi  +  Oi< 

LUO< 

mm-*-jN*-<*- 

<x 

1-4 

a>— «- 

-^ 

■~<M-0-* 

— . 

O— >« 

—  X—  s 

o 

»  _   <u  _i  t 

t— 

ti- 

•■* — 

U- 

-*>>> 

u_ 

•— * 

-JLUUJI—  — ' 

ujoc^o 

OOOOOOOOoOOO 

UJU. 

ll. 

C\JO  O    1 

u, 

>       * 

U_ 

cMrgO  | 

<-#23rs# 

acuji-i 

X-. 

tn^>*- 

UJ 

w  Q> 

UJ 

*  +  +fM 

UJ 

w>0 

x—  t—  Oct 

XX 

1  1  1  1  1  1 

H-U- 

OOOOOO 

a 

+    •  +  — 

o 

o         < 

O 

+  +    .-. 

*  -JO*  Lucom 

U.UH 

C\JU\U> 

o 

r\ic\i>rm 

u 

f*1    .    ■»      ^W 

o 

vU-^fMU^ 

x«-aoi-uzo 

o~ 

U-LU 

OOOOOO 

.>^.;>;> 

•  f\|>U 

3»>^*> 

2TX  1     •         »CC<X. 

Xai 

cm  cm  mm  mm 

ox 

v»v.*n» 

—4 

www  i 

CNJ 

CM>>   I 

CM 

WW   | 

>— «— wm  ii  -*uj 

oo    a 

h- 

co-^**" 

• 

*  **-• 

• 

*** 

• 

**■*-* 

H.  II 

UJ 

UUOU. 

oooooooooooo 

00 

OOOOOO 

m 

-^.-4.    >> 

ro 

CMCMfM-^ 

m 

mmm> 

II  II  II  II  —  II 

D 

uz 

OLU 

(\jmm 

<<<t^- 

<r«> 

<<<^, 

-5       II  — 

Z 

Z<UJ 

il  il  il  H  ii  il 

— •»— 

OOOOOO 

z 

wwwv 

Z 

ifUIUH 

z 

^»^—  >M»«^ 

ct:LU-*<N— ro      -j 

»— * 

1X1     s 

l-< 

a 

o 

o 

araczsrssrs-— w 

— 

acuj<t 

coin  — *~4*t  -*■ 

<z 

II  11  11 

H-4 

II     II     II     II 

— • 

II  II  II  II 

►— 1 

II      It       H      II 

i—  i—  2  oc  oa:  ->o  «£ 

lUCCOO 

OOOOOOOOoOOO 

QCf-l 

1— 

H- 

f- 

ooiuujziU'-za 

U-<1 

cMCMcnmcnm 

s: 

~4c\jcn 

< 

— <c\i  en— < 

«a 

-TinsOfM 

< 

i*-cocj*m 

oat—  t—  <h- oocto 

LL3UJ 

oo  on  oo  oo  on  oo 

UJf- 1 

<<< 

o 

ooo**£ 

O 

ouo>: 

o 

OOU>i 

"UI 

J- -J 

o 

car 

CJ 

o 

QOO»- 

V-UJ 

UJ 

UJ 

UJ 

<\J 

uuuoo 


uoou 


uoo 


ooo 


ooo 


46 


oooooooooooooooooooooooooooooooooooooooooooooooo 

>f  <4"  «4"  «*■  -4-  it*  in  in  in  in  in  in  in  in  in  o  *o  *o  o  <o  vo  ^  <o  -o  o  r-  r-  i**  r- 1**  f*"  r-  (*-  r»»  f-  oo  oo  ao  ao  oo  oo  ao  co  co  cncroo^ 

oooooooooooooooooooooooooooooooooooooooooooooooo 
ooooocoooooooooooooocooooooooooooooococooooooooo 

UJ 

u 

h-  Z 

2:  < 

uj  — 

CC  00 

UJ  — 

U_  Q 

u_ 
a 

<  z 

o 

u.  •-« 

-•  I- 

•  00 

2  a 

a  o. 

►-  •  x 

>—  LU  - 

OOCL  •• 

QUI 

za 


• 

INI 

UJ 

INI 

z 

«» 

>-* 

> 

K- 

> 

X> 

•» 

O 

X 

ac 

X 

•JQ 

** 

3 

m 

OO 

X. 

UJ 

c* 

_J 

0 

13 

•» 

oc 

ao 

O 

00 

•» 

• 

h- 

Of 

O 

UJ 

•ft 

3: 

CVJ 

< 

•*. 

cC 

m 

0 

vO 

0 

UJ 

•ft 

X 

in 

<o 

UJ 

0> 

OO 

_J 

0 

< 

H 

u_ 

• 

O 

•» 

OO 

m 

• 

•> 

o~ 

<X 

ox 

O 

MU 

O 

o 

Q      -*0  ao  <m  ao  O 

►—     Oao  ao  I— 

a  o 

h-    — <o       >—  a       f—  o 

uj     «.  a  00 

00    mo       000  •  — 

UJ        *  —  -»  O  Q 

QC      CJ—  O—OO  • 

<    00  •  o  •  in  o 

•*  •        o  •        o  o     <*-•  *—  a 

OO      -10  OOO+  O  00  H-  — *— »-.—.— »—  t  X  cD 

h-  00        000                            m  oouj  +  o  o^^^i^:^^:      •-»              a 

z  ■— o        ^o-*     o        mo- •     — »— »              •  >— iq  _j  o— ■—»—«—'——— >—  i£     x        o 

uj  t—ino      1    •  '-«     mo      1    •-*     fH^         2:uj  -»  u»  (nx>ccx>n^j-     o     co— 

—  <    •    •        »0*-X    •    •        •O*- >"—  <—  INI          CD  Si/)  — •  oCtCCQiUOO—    •>      o      ^Q 

O  3!t— O      hOXXKO      H-0>>inIMin|-*  II     •  aC—»  —  UJ  -4  I    |     |    |     |     |  >-»0      X      ^Ol' 

•-*  OQjO-^OO           JOWO                                 —  O  niO  X  II                               I  n        p-      uj^iO 

a.  z>  •oco  «ir\  11  11    »ooo  «in  11  11  11  11  11  +  —  -«  u_a.         m  >—  >£  11  11  11  11  11  11      -»  11  cr  n  ai 

u.  ooxm     x  1         >m     >  i                       2—  h-uj  h                                         11  >            zi— Q- 

uj  x    ax    —- •>    c>    -.—-»-».-. Zwivi  <f-        a  00  r--.— ^-. .    ^~.ci~.*+ 

Q      jw  ||  |-w  ||  ^Zw  11  |-w  11  2ZZ221      t— ^-      —JO      —  •—      ~      C^^T^^iii^^X— ^1— vf»— _i-J 

O        -J  «— «-  ^.^..-^.w  ||  M  Q.  — (  — w  — —  0        _        wZ_|_J 

<U.XOU-XXXLL>OU.»>NQCN       UJU.  UJ       wO       a)       OX>CiX>NOlLOOUU« 

uj    o^xo^xccoi-ixjj^-xaiOoccci-jziNii-i  ujini  no    hh  Qaccocoooi'-'>o>uoo 

a.  j_i-i  x 

h-                   o     *-*               cnj     m  1— 1/>  h-  o                                     or- 

00     ao  00     00  >f  cj^u^ 

uuu  0000        000 


a.  00 

uJ3 

1—  *■> 

oOQ 

< 

UJ 

Suj 

1— <TQ 

^« 

►- 

• 

UJoO 

UJ 

ZD 

Z 

az 

1— 

ai- 

x> 

UJom 

a 

O 

O 

oi 

z  » 

a 

<Q 

• 

CJ3 

>UJ 

0 

z 

Qa: 

m 

>— » 

<-* 

*— 

to 

H- 

OOUJ 

+ 

O 

—•a 

— 1 

-^ 

u. 

2:00 

— ( 

CC  — 

— 

UJ 

O 

NlO 

X 

u.a 

m 

1— 

H-UJ 

11 

<h- 

O 

00 

_|CO 

-^t— 

— < 

Q. 

— H 

UJ 

wO 

00 

UJINI 

NO 

»-« 

47 


oooooooooooooooooooooooooooooooooooooooooooooooo 
m  *-m  vO  h"  oo  0*0  — icMm  >$•  m  >o  i**  oo  o  o-^  cMco^-m-or*-  oooo^r\jfO«*ifv  *0f**  oo  o^o-^rMm**^  sOf"»oo  o*o 
O^o^^^o^OOOoOoOOOO— «.-^^^^-*^^— rf^tM'NcgfMCMCMfNCMrMtMc^rnrOr^ 
— *  .-»  -4  -^  ^  -<  *-•  eg  r\i  (M  cm  c\j  cm  rvi  eg  cm  cm  cm  cm  cm  cm  eg  cm  cm  cm  cm  cm  cm  eg  rvi  cm  (\j  cm  eg  eg  eg  eg  rg  eg  eg  cm  eg  cm  eg  eg  eg  eg  eg 

oooooooooooooooooooooooooooooooooooooooocooooooo 
oooooooooooooooooooooooooooooooooooooooooooooooo 

••f—      LU->—UJ-l—       ->*-  OO 

•        000-LOO-00****  << 

«A     •— *      "Z      *&     — •      Z     **     •— i     Z       "     Z  at  QC 

Za<ZO<ZaOa:0  UJUJ 

0->—Q«»»—      0"1— <CC»-»  >> 

_4»00*--00t-t-H-OI—  «^< 

I—        -  .— i        ^—        m  — .        H-         »  hh        (J        « 

i-*<I*>Qi-*<AQi~«<AOO-</)  LUUJ 

oOZ-lOZ-i/>ZO-0  XX 

00*00*000.-0.  H-l— 
Q.      ■— •      •        Q.      •— i      •        Q.      •— t               «*> 

h-      •**>  I —      *»^->Q>  oooo 

XwZ>hZNh  UJ  _,-,.-, 

.     ooo-Loa-Loa»—»  - 

»Om          »Qm         »OUl<(i  -  -  • 

-  Q.  t—  -  CL  h  ■  Q.  h-  -J  —  GC  ■»  •»  LU 
**  — .  &  m  4A  O  3-  O  '-•  »-l  3 
OXoOO>oOQr>giU      0**0  •  •  -J 

lu     -      o     lu     •      o     lu     -      at     JA»  p*        r«-        < 

H--0.1—      -o.h-      -at<i-  u_u_> 

<t»                <      -                <      -       U     V-M  •  »  •> 

_J**>— J<*r«g.J«*OI**  X  X  ivl 

O— J-O— l»3—l-»— J  — «  — t  — 

U<       »U<       »U<       *     OCX  >v*n,Q- 

-J     O     -       -JO-       -JO-       at-  O  lu  • 

<►-**<(—      *»<►-«*      O   H—  •  -  H-  O 

o      o               O      <->               O      O          —O-  O***  LU  LU  <  Z 

-  -»<     •      •  — <     •      •  -»<     .  «  »  «»<r»  o        o        -J    —« 

•>•  -  ••        »•  ■  »        mb  m  -**    -A-   *»  Z  Z  O      O 

.    -*    »      -        .  <a   •>      •        -  •**   •>     .  A-     I     -A  <  <  O      < 

ZA»        Z      ZA.        3      ZA.        ZIZI-    I  I—  H-  -I      LU 

M|Z  M  «|Z  -  M|Z  H|M|Z|  00  OO  <  X 

-I     I     •"•—-J  -I     I     •-—-J  — )     |     •——-J     |     _J  HI     |  _|  _  O 

Z  I  _J-  Z      Zl   -J»Z      ZU»Z      Z*_J  O          O  51 

— «       Z**'-'       i—»       Z***-1       h_       2<*N'>wOZ>  LU  C£ 

-JX**A_|      — l>~>A_j      jKmAj»  _!<.-..  x           >          o  a 

-»«  -  _j  i  •  —  •  •  _j  i  •  -*•  -  —j  I  -    -•  at— i  -  <  a. 

■      «   »•     |      kb      •    ~»     |      »»      •.   •—     |      •.»      «.K-»    -  LU  LU  QC  h" 

**    •■      •>  |      «*A    •  -      ••  j      »4A    ••      »  j      **A    •         »-^  O  O  LU  ^ 

ALU**    •       LUAUJ«A    •       lUALU<*    •       JJALUO    »A  Ot  at  >  -J 

I   OOALUXOO   I  O0ALU>00  I   OOALUMOO   |  OOZLU   I  F— -»  f— -»  <  Q. 

i  _i  i  o-  j  i  j  n*  j  u  I  o-  _i  i  jo-)  i 

I  <  i  oc  -<  I  <  I  a:  «<  Klan  -<  I  <_ja£  I  lu  •     lu  •     lu     lu 

LL  I  >--  LL      U_  I  t—  -  U.      U.  |  t— -  LU       LL<>-  XLU       XLU  X  X 

X    •         •<*>    •>•    •         •**    ,,SJ    •         •■**    »^    •         »^  f—  O      *~"0  >—  K— 

•     -^  *A   ••-     »^:   -A   •»-     »^   «>A   -O  ••LU  m_)  2          z 

-ooo  I  m  -ooo  i  m  -ooo  I  m<ooo<  xooooo<  oo  oo 

-  o<o  i    —  o<o  i    —  o<o  i    »oiozoa:  o«h-o>— i^-  •-•  »-• 

*»rga;rM  I  <**rga:cM  I  <^<M(Xcm  I  <H-(M<fM>—  >     oo;>     oo 

A    •"♦—   -      OA   -H-   -      OA    •+■■   -      O        -t—    -  «-•-»—.  «t.  -<<.  — jvj-  Q« 

I  x     x^o  i  >    >*:o  i  rvi    oi^oooooxj)  ^^^  -  -a  *  -oo  -<i  - 

I  OOQCO   *  I  UOQCU    -  I  OOOCO   -Z>— .a:^:  O—— —  X->      >-^      >-4LU— • 

I     -2   -<CQ  I     -Z   ».<CO  I     -Z   »<COO   -O    -O  OX>M  OOO                I—    «X»—    •><    »X    • 

rvJO^QCO     NQNoia      ^OivlQiO— !•*>•"  ^-J  (\jOUO  000                   p^          r-          r*»      r- 

^.cc  —icth- Oi^ac— iqch-o^cc-jqci—  o<icx  -ot;<i  •>  ogrMfM^         — luo— u.a— u.—u_ 

0«-»<l*~     «-o>— «5^«.     ,w ^  — <x>—     ««     «»■  ^*      _j  .^^  ^  >  wno          tr\   -ujo  -ujro  -T*-  - 

<Q      QOQ<Q      000<0      QOQQJQ«Qujii  >>>vO—»-»Q0XI—  COXJ-r-XCOX 

QO-UJHZKQtll-UJt-Zt-Qtt-UJKZH-UI-ZH-OUJOOO  II  >»  <<<    --*'•*   •— )<    •— <<    »— ;   »^ 

i-auoooi-ooOQOhQuoaozoQozz  •  •  •    <<<Lux>tvj^*——  ^o— '-jvo«-_jvo«->o^- 

JZ- 1— l-J      — IZ— I I  — I      _lZ-l-)-/<-l'-_«OOOO^X>ND  «-X>-'-.      o^-      o^»     ~» 

00.<0.<0.00.<Q.<Q.OO-<0.<CLf—O.H'0.l— Q  Z  II    II    II         l|        HO      I— O      ►-      H- 

Z       I—  Z       h-  Z       H  00       m       OO         II    II    t|  m  II    II    II  —<  LU  LU<_JLU<_JLU<LU< 

0_ioo-jlu_jo_joo-jlu_jo— ioo_ilu_ji-._joo_ih-,_j  ^  >—  x><^if—  11  11 1— 5-<r^y<H-T-H-s: 

— !_)•-•— ilj_j_j  _jt-«_jo -J  _j— '»-•— jo— jo— iq_iq_i;>;>>    >>>zooo«-i        m •cco'—QCO'-iQt'-'Q: 

«a<Z<«o<Z<;i<i<Q<z<»  <o.<»  <<t<<tr_j<<«io>^>>Qix>a:o     Qto     ctOa:o 

O-  LJ<0      O"  0<0      O-  i-)***^    -O      O    •OX>NQX>NO«<^I-I-3U.UJ3:U.UJ3U.3U. 

m  m        so        rn    r- 


48 


M 

rsl 

CO 

CD 

CO 

CO 

* 

•* 

# 

> 

< 

> 

< 

< 

< 

~^ 

— 

-^ 

— . 

>— 

■>^ 

* 

# 

* 

O 

X 

X 

o 

o 

o 

ooooooooooooooooooooooooooooooooo 
— •  rg  m  >f  m  so  h»  co  o  o  «-•  rg  rn  -^  in  «o  h»  oo  o>  o  -^  eg  m  ^  m  o  r»-  oo  a*  o  — •  <n  m 
<•$■  «^  >*•  -J*  <f  »*■  «-f  -4"  <j"  in  lti  in  m  in  in  in  in  m  in  <o  -0  >o  ^o  *o  <o  >o  >o  ^o  vO  r»  h»  r-  is. 
rg  rg  rg  rg  (\i  f\j  rg  <\j  rg  cm  nj  eg  t\j  rg  (\i  eg  (\j  rg  rg  rg  cm  oo  rg  c\j  rg  f\J  <\j  rg  rg  rg  <N  rg  rg 

ooooooooooooooooooooooooooooooooo 
ooooooooooooooooooooooooooooooooo 

M  <      <      < 

<  w       w       — 

—  Ill 

rvl  -» 

ivl  f>i 

»  CO 

>-  # 

>  > 
< 

X  — 

X  — 

•o    # 

UQ      X 

»^>g 

rvlfvl       + 
O    » 

»x  »o 

XXO-H- 

O  -Ujrvl 
•Oh< 
COO<W 
CO    »_J  I 

CO    .JU 

CO    •»       < 
•X'-U— 

XUCQ— 

2T  CO    »       ■* 
<00  •COV-X 

»-•  <tcQOcO 

Ql  <  -Z  III 

l-4«  «Nl2    I 

rvlcO<  -»-»—. 

<    »0—  —      —      «■« 

»>      —      00  >      o      > 

>cox>    0^        o    o    u 
uj    <  •••-•o  *     #     * 

Z        »Xa£#       O  fsl      r»g      CO 

m*       XcOh-rvj       h-  CD      CO       CO 

H-         <     »<CO  —         —         w 

—2:        -   >s  o    —<s^    a       111 

a.-  -  a:-»os.  -  a     1— <•   1      O        «■»-»-• 

•— •   •>    >00    •q.    •.  •  c£      <T    ..  •— »  M      rvj      o 

Q-^ZI    •  O'— >-*Z  CO      ST^JXrvl      ~»  <*J      >*J      <*J 

uj  «:s    OvOQ  »^  o     co<vOO     o        *■;(■■* 

r*0  +  v0m     r>-C3  1/)     o  ••oj*.      Q         >     co     > 

-»U»C_3      rn      <^U_Q      -^.<—  tO>    »>.         •  CD       CO      CQ 

CO    ••      U    •   I  U*    »  O*  LU  <I  m   CO       O—       **      ""^      ■«* 

COXl/)    •  »—       COXl/>       -O  _J       LU    »  +  «—  »r»-       *-Q  —  Q  — Q 

vO—     sO  •'.UvOw         so—  a:     ►-•<—*      Z^0O<nxnxn 

*-     +m(ja*.     4.     w  (_        x       »^.o^<— <t— <  — 

t—  LU  I—  t—  00      3C0H-<      Q  w— w^»^,-»2 

LU<    *IIQIlLU<    •»      LU<  •         O*  <       — 'QUJO       —      -«       — .CC 

^S0.  ."-  .'"St***-^11-^     <£     o;_jsr  u  — — »-i— 11  v  11  co  11  >^ 
•—qc-jo     OwaZ»  -*aczDd     lu     co<o:     >-     —•         ca     co     cq>-q 

QCOLULULLLUQiOLU  «o:Oh-k:  SI  -3LUOQ(DU.CtOX*>#N*lU2 
-SU_»—  Q«Q3U.SlilSu.olLU      <      i^OQiXQ* '—30Xrv4>'v|fv)<tCXLU 

*  **  at  *  *■** 

CO  ^JTJ*  OO^  CO  l/» 

co  coco  g3  0s  o^ 


u_ 

OS 

SK. 

0 

<S) 

H-» 

UJ 

H- 

-JUJ 

LU 

ax 

z 

0.1— 

0 

•-* 

< 

QU. 

3: 

O 

00 

«>o 

•  1— 

• 

•2T 

•» 

•  LU 

• 

IZ 

X 

-40 

«-4 

*a. 

•i 

-  2: 

■ 

t-a 

X 

LUO 

h- 

C3 

ar 

a:  _i 

< 

<< 

UJ 

i— h- 

z 

LU 

LUO 

X 

IN 

1— 

1— >-• 

ac 

U. 

LUO 

a 

OX 

LU 

LUUJ 

-J 

-JX 

CD 

CJ>I — 

>- 

O 

> 

<-) 

O 

0 

* 

* 

# 

rvl 

rvj 

< 

<5 

< 

< 

I 

INI 

^g 

0 

O 

O 

0 

* 

* 

-* 

> 

< 

> 

< 

< 

< 

^* 

** 

^» 

•«• 

«— 

•M 

* 

■!<■ 

# 

CO 

X 

X 

CO 

CO 

CD 

O 

^4- 

1-5  * 

CO 

— « 

LU 

LU    • 

X 

O 

xr*- 

h- 

I— 

H-U. 

0 

00 

O 

LUl/)X 

•— 1 

L3 

Qw^ 

49 


APPENDIX  B 


COMPUTER  PROGRAM  FOR  NUMERICAL  OUTPUT 
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APPENDIX    C 
EXAMPLE    OUTPUT    OF    APPENDIX    B 

An   example    of   the  output   of   the    program   in    Appendix   3    is 
shown   below.  The    total   output    for      a    single   time      st ep   is 

presented. 

SENSOR  SLANT  SENSED  CALCULATED 

NUMBER  RANGE  SIGNAL                           SIGNAL 

5C50.Q  0,60540-04  0.6054c-04 

2  506C.1  0.1596D-C3  0.1612E-03 

3  5050,2  0.26820-C4  0.2683E-G4 
A  5099.0  0.3940D-C*  0.4135E-04 
5  5050.2  0.9366D-C4  Q.93b8E-04 

ECLATIGN  TRUE  APPROXIMATE 

S2-S3  0.1327C-03  0.i.34j>E-03 

^2-S5  0.6590C-04  0.t>7*8E-04 

S3-S1  -.3372C-G4  -.3371E-04 

S5-S1  0.3312C-04  0.j314c-G4 

c3-«;4  -.1258004  -.it52E-04 

S5-S4  0.S426C-Q4  Q.5233E-04 

EQNS                                                                   X  V  I 
USED                                                                 POS                   POS                    PCS 
12    3            ACTUAL                                      0.0         -1000.00       -4950.00 
C^l^LAI^a z&a*.£.8. z£LZ*&L—z<tli<l^l 


The  first  fi/e  rows  of  data  above  are  for  the  appropriate 
"sensor11  as  indicated,  the  slant  range  from  that  sensor  to 
tha  target,  the  received  signal  as  generated  by  equation  2.5 
and  the  signal  that  wculi  be  generated  using  the  approxima- 
tion described  in  Chapter  3  to  linearize  the  problem.  The 
next  six  rows  show  the  results  of  taking  the  differences 
between  the  received  and  calculated  (for  comparison  purposes 
only)  signals,  listed  under  the  headings  of  TRUE  and 
APPROXIMATE,  respectively.  The  final  two  lines  shew  the 
actual  coordinates  of  the  platform  relative  to  the  target 
(the  coordinate  transformation  has  not  bean  done  at  this 
point)  as  compared  to  the  calculated  position.  These  are  the 
points  used  by  the  program  in  Appendix  A  to  plot  the  graphs 
after   the   coordinate    transformation    is    completed. 
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The   second    section  of  the   output    is    shown      below 


EQUATION 

CJ/SiSTANT 

X 

NJMbcR 

COEFFICIENT 

1 

G.5  JD+G3 

-0.30Q+G2 

2 

-0.33u+03 

G.  14u+GG 

3 

-O.SoG+02 

-0,210+01 

Y  Z 

COEFFICIENT    COeFFiCIeNT 

-0.920*01  G.2GG+0i 

G.  130+QG  u.  54J-Ux 

-0.190  +  02  G.33lh-01 


^ATIO 


1/2 

-0.13E+G1 

-0. 22E+G3 

-0.  71c* 02 

0«  3oc*-02 

1/3 

-0. 3 6c* Ox 

0.14c>G2 

Q.t9e*00 

0  .56c+0Q 

2/3 

0.  44E+G1 

-0»  c5c-0x 

-G#  o9c-02 

0*  15c— \j! 

The  first   three  lines  of   the  output  represent   the  ccef 
cients  of   the   equations   3.6   -   3.8   in   the   form 
CONSTANT  =  X  ZOEFF(x)  +  Y  COEFF(y)  *  Z  COEFF(z).  In  order 
determine  the   linear  independence  of   these   aquations, 
ratios  of  the  coefficients  are  taken  in  pairs.   Since  it 
be  seen   from  the  last  three  lines   of  the  output   that 
ratios  of  the   constants  and   coefficients  differ   signi 
cantly  for  all  three  equations,  it  :ai  be  concluded  that 
equations   are    linearly   independent.    This   holds   t 
throughout  any  given  run,   which  may  be  verified  by  print 
out  all  time  steos  for  that  run. 


ri- 


cf 
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ing 
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APPENDIX    D 
CRAMER'S    ROLE 


A.       6ENEBAL 

Cramer's  Rule,  a  method  for  solving  systems  of  linear 
equations  in  two  or  more  unknowns,  states  that  each  unknown 
can  be  expressed  as  the  ratio  of  two  determinates.  The 
general   form  of    a  solution    is    shown    below. 


Ax+By+Cz=0 

Ex+Fy+Gz=H 

Ix+«Jy+Kz=L 

ABC 

D   = 

E    F    G 

I   J    K 

0    B    C 

x    = 

H    F    G 

L    J    K 

A    0    C 

7    = 

E    H    G 

ILK 

ABC 

z    = 

E   F    H 

I    J    L 

/   D 


/   D 


/   D 


B.       EXAMPLE 

The    following  example   v as   taken    from   [Ref.    4]. 
3x-2y+2z=7 
x*    y  +    z=6 
2x-    y-2z=2 


Solve 
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Let    D   be  the  denominator 


D 

= 

I     3   -2      2     ! 
i     2   -1   -2     | 

s 

-17 

X 

= 

|    7    -2      2     | 
I     6       1       1     | 

LLrl-iLl 

D 

= 

-51 
-17 

y 

= 

I     3      7      2     I 
|16       1| 

LL.LzLi 

D 

= 

-34 

=    2 

-17 

z 

s 

13-2      7     | 

1       1      6     I 

12-1      2     | 

= 

-17 
-17 

CHECK: 

3x-2y+2z=7 
3(3)-2(2)+2<1)  =7 

9-4+2=7 
7  =  7 

x+y+z=6 
3+2+1+6 
6  =  6 

2x-y-2z=2 
2(3)-2-2(1)  =  2 

5-2-2=2 
2=2 

C.       SUBROUTINE    "CHAMBERS    RULE"    CHECK 

Running  the  values  in  Section  3  through  the  subroutines 
of  the  programs  in  Appendix  k  and  Appendix  3  yielded  the 
same    results  as    the    examole. 
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